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Q: Abstract 

Da In the next-to-minimal supersymmetric model (NMSSM) a light CP-odd Higgs boson is so far 
= allowed by current experiments, which, together with a large tan 8, may greatly enhance the rare 


dileptonic decays B > X,0*(^ and B, — (* (^^. We examine these decays paying special attention 


" 


to the new operator allowed by the light CP-odd Higgs boson. We find that in the parameter space 


Me «n 
"Hu —E 
e: 


allowed by current experiments like LEP II and b — sy, the branching ratios of these rare decays 


can be greatly enhanced and thus the existing experimental data on B => X;u*pu- can further 
stringently constrain the parameter space (especially the region with a super-light CP-odd Higgs 
boson and large tan 3). In the surviving parameter space we give the predictions for other dileptonic 


decay branching ratios and also show the results for the forward-backward asymmetry. 


PACS numbers: 14.80.Cp, 13.85.Qk,12.60.Jv 


I. INTRODUCTION 


Recently some non-minimal supersymmetric models such as the next-to-minimal super- 
symmetric model (NMSSM) have attracted much attention [1] since these models can solve 
the p-problem and alleviate the little hierarchy. In the NMSSM, for example, the p-term 
in the superpotential is forbidden by imposing a discrete Za symmetry and instead it is 
generated through the coupling between the two Higgs doublets and a newly introduced 
gauge singlet scalar which develops a vacuum expectation value of the order of the SUSY 
breaking scale. In this way, the u parameter at the weak scale can be naturally explained. 
The NMSSM can ameliorate the little hierarchy by either tuning the parameters to enhance 
the theoretical upper bound for the mass of the lightest CP-even Higgs boson or relaxing 
the LEP2 bound of 114 GeV through allowing for a light CP-odd Higgs boson (41) with 
mass below 2m, |2].. 

It is interesting to note that in the NMSSM the lightness of such a CP-odd Higgs boson 
can be naturally predicted in the enlarged parameter space, and is also allowed by the LEP 
II data li. This light Higgs boson can not only alleviate the little hierarchy, but also can 
help to explain the observed anomaly in the decay X* — pu*pu” |3]. On the other hand, if 
a Higgs boson is indeed so light, its effects in some low energy processes may be sizable and 

ü : For a super light A, the decay b — A,s is open and an 
4 


thus are necessary to check 
(note that as analysed in lh, a small CP-odd Higgs mass 


analysis has been performed in 
is only protected from RGE-effects in the limit of large tan 8). In this work we consider the 
full possible mass range of A; (heavy, intermediately heavy and light) and check the NMSSM 
effects in the rare B meson dileptonic decays B — X,/*(- and B, — (*€^»y |6j. 

These rare dileptonic decays are induced by the flavor-changing neutral-current (FCNC) 
b — s transition and are of special interest due to their relative cleanness and high sensitivity 
to new physics. In the Standard Model (SM) such FCNC processes are suppressed and have 
very small branching ratios but can be greatly enhanced in some new physics models 
ld y id Al hd id. Since experimental data on B > X,u^*p- is available and the future 
LHCb or super B factory will further scrutinize B meson decays, these dileptonic decays 
serve as a good probe of new physics. 

In supersymmetric models these dileptonic decays can be drastically enhanced by large 


tan 8 since both the b — s transition loops (such as the charged Higgsino loops) and the 


Higgs couplings in the Higgs-propagated diagrams are proportional to tan P. It has been 
shown (see e.g. |12]) that in the minimal supersymmetric model (MSSM) great enhance- 
ments are possible for these decays. In the context of NMSSM, in addition to the tan 8 
enhancement, the presence of a light CP-odd Higgs boson could further enhance these dilep- 
tonic decays. Thus, the parameter space, especially the region with a super light CP-odd 
Higgs boson and a very large tan B, is constrained by the existing data on B — X,u*p-. 
In our analysis we will examine the NMSSM effects in these dileptonic decays by scanning 
over the parameter space allowed by the LEPII experiments and the data on b — sy. We 
will show the 20 constraints from B — X,;u*p- on the parameter space and then give the 
predictions for other dileptonic decay branching ratios and forward-backward asymmetry. 

A key point in our calculations is the presence of a new operator due to the light CP-odd 
Higgs boson. In contrast to the MSSM, where all Higgs bosons and sparticles are so heavy 
that they can be integrated out at the weak scale, the light CP-odd Higgs boson A, in the 
NMSSM cannot be integrated out at the weak scale and thus a new operator O4 describing 
the interaction A,bs must be treated carefully. 

This work is organized as follows. In Sec. a brief description of the NMSSM is 
presented. In Sec. MI] we calculate the Wilson coefficients paying special attention to the 
new operator O4. In Sec. we evaluate the NMSSM effects on the dileptonic decay 
branching ratios and the forward-backward (FB) asymmetry, and present some numerical 
results. The conclusion is given in Sec. |V| and analytical expressions from our calculations 


are presented in the Appendix. 


II. A BRIEF DESCRIPTION OF NMSSM 


In the NMSSM a singlet Higgs superfield $ is introduced. A discrete Z4 symmetry is 
imposed and thus only the cubic and trilinear terms are allowed in the superpotential. The 


Higgs terms in the superpotential are then given by 
ASH Hp se. (1) 


Note that there is no explicit jj-term and an effective -parameter is generated when the 


scalar component (S) of $ develops a vacuum expectation value s/ 2: 


Heff = A(S) = (2) 


um 
pol] 


The corresponding soft SUSY-breaking terms are given by 
Hie. Els 
A ASH, $ Ay + q i (3) 


The scalar Higgs potential is then given by 


Ve = |MH,- Ha + KSP? + JASI? (|Hal? + | Hs) , (4) 
2 2 
g G 2 
Vp = > (|Hal?|Hul? — | Hu - Hal?) + = (Ha? — |H, 15", (5) 
Vos = màéHa[ + m2|H,P + mls? + (ASH, dra AnS + h«c.) , 6) 


where G? = g? + g2 with gı and g> being the coupling constants of Uy(1) and SU;(2), 
respectively. 


The scalar fields are expanded as follows: 


1 . 
-= (vat da tt Hj 
mel RT D yy. — da 


l (s+0o +i). (7) 
Hz Jg (Vu + bu + ipu) 


Sl 


The mass eigenstates can be obtained by unitary rotations 


HA; Qa 4 li Yad 
H Gr Hy 
Hə U Qu E A» U Pu E U ) (8) 
H* H* 
Ha oO Go £ 


where H123 and Aj» are respectively the CP-even and CP-odd neutral Higgs bosons, G^? 
and Gt are Goldstone bosons, and H* is the charged Higgs boson. It is clear that the 
charged Higgs sector is the same as in the MSSM, while the neutral Higgs sector contains 


one more CP-even and one more CP-odd Higgs boson. U^ and U^ are given by 


Co,Se CoCo Se, rang (Con — 20404) Con —Son 
U^ = — So, SB — So, C3 Co, > p = tang (Son + 264Co,,) Son Coy 2 (9) 
-Ce Sg 0 1 map senpi 


where Cy = cos X and Sx = sin X (X = 04,05). The mixing angles are given by 


T U 1 2Aus 
Se » E 2 = 1 
Ya 2 T ad n (=) y an 24252 — M? eu 
with v ~ 246 GeV and 
2 
y - PARIS (11) 
V/2.45 + KS 
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The Lagrangian of the Higgs couplings to quarks for a large tan 9 are given by 


E = a (56 Ai, tan 342) dysd, (12) 
i 2mw S 
gom, 1 Ó v 7 
L ¡Qu — — ———_A a > 13 
m "2mw tan 8 iu : ) e us) 
goma U U I 
Lun mmy = 55 ((Co — 2848) Hi, (So + 26.C9) Ho, tan BH) dd, — (14) 
My S S 
gom, H; NY. 
Lies 2mw (com. SoHo, y) uu. (15) 


Here one can see that the coupling of the neutral CP-odd Higgs A; with up-type quarks are 
suppressed by a large tan 8 and thus can be neglected in the large tan 9 limit. 

Since one more Higgs superfield S is introduced in the NMSSM, we have a new neutral 
Higgsino wg. So the neutralino sector is composed of Uy (1) gaugino A!, SU;(2) gaugino X, 
and the Higgsinos Vi, Va, and vs. The corresponding mass terms are given by 

Lng = EUC? — gA!) (cos Busy, — sin OY) — SMN — EMAA! — Asi Vd 
t 2 d u 2 2 vo a" Hu 


GA (cos Bt Us + sin Bibl bs) + yen + h.c. 


= SQ Y vw + he, (16) 


where 


(997 = (—iA!, iA, Wh, Vir, Ws): (17) 


The neutralinos are obtained by the unitary rotation y? = (Z Ni, where Zy diagonalizes 
the mass matrix Y\o. 

Similar to the charged Higgs sector, the chargino sector of the NMSSM is the same as in 
the MSSM with y replaced by ness. The chargino masses are obtained by the diagonalization 


of the mass matrix with two unitary matrices Z_ and Zy: 


-M> V2my sin 3 
v2mwcosÜ —Hesf 


III. CALCULATIONS OF WILSON COEFFICIENTS 


In our calculations we consider the flavor mixing between b and 3, which make contribu- 


tions to the dileptonic B meson decays through gluino or neutralino loops. Following the 
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analysis in , we assume the flavors are diagonal at tree level and the mixings are induced 
at loop level. Such mixings can be parameterized by a small mixing parameter e, which is 
dependent on some soft-breaking mass parameters |15]. In our numerical calculations we 
input e, = 0.1 for illustration. We perform the calculations in the Feynman gauge and thus 
the Goldstone bosons will be involved in the loop diagrams. 

Since in the NMSSM the lighter mass eigenstate, the CP-odd neutral Higgs boson Aj, 
can be rather light, with a mass ranging from 100 MeV to the weak scale ill ( 
package NMHDECAY , we checked that such a light CP-odd Higgs boson A, is indeed 


using the 
allowed by the LEP II data), in our calculations we pay special attention to this wide mass 
range of A, and discriminate three cases. 


(i) Case A: Heavy Aj. 


For a heavy A1, around weak scale, we integrate it out together with the other heavy 


particles (Higgs bosons, top quark, W= and Z bosons, sparticles) at weak scale to ob- 


tain the Wilson coefficients. The effective Hamiltonian describing b > sI*I7 transition 


reads 
AG 10 2 
F * 
Hay = =- 7g VeVi 2, Cilh) Oil) + D ,Cou)Q(u) . (19) 
i=1 i=1 


where O; and Q; are operators listed in d "n and C; and Cg, are respectively 
their Wilson coefficients, and ju, is the renormalization scale. Note that the most 
general Hamiltonian in low-energy supersymmetry also contains the operators O! and 
Q' which respectively are the flipped chirality partners of O; and Q;. However, they 
give negligible contributions and thus are not considered in the final discussion of 


physical quantities L 


In this case, there are no new operators and we only need to calculate the NMSSM 
contributions to these coefficients C; and Cg, at the scale of my. For the processes in 
our analysis, only C7,9,19 and Cg,, are relevant. The NMSSM contributions to C7 9,10 
are the same as in the MSSM, which are computed in 18 . For CQ, , the NMSSM 
contributions are different from the MSSM contributions and thus need to be 
calculated here. The Feynman diagrams we need to calculate are shown in Fig.1, 


where the loops respectively involve the charged Higgs bosons, charginos, gluinos and 


FIG. 


neutralinos. From the calculations of these diagrams we obtain the Wilson coefficients 


at mw scale, which are presented in the Appendix. 
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(h) (i) (3) (k) 


1: The Feynman diagrams which give the dominate contributions to CQ, ,: (a-c) charged 


Higgs loops, (d-g) chargino loops, (h-k) gluino and neutralino loops. 


(ii) 


For the calculation of the dileptonic B meson decays we need to know the Wilson 
coefficients at the m, scale, which can be obtained from the running of the coefficients 
at the my scale down to the mj scale. Such a running is governed by the anomalous 


dimension which can be found in ] 


Case B: A, with an intermediate mass my « Ma, < my. 


At the my scale we retain A; as an active field and thus we have a new operator 
O4. After integrating out all heavy particles at the mw scale we obtain the Wilson 
coefficients including C4 for O4. Then we work out the anomalous dimensions and 
run the Wilson coefficients from the my scale to the ma, scale. At the ma, scale we 
integrate out Ai, which makes an additional contribution to Cg,(ma,). Finally, we 


run all the Wilson coefficients from the ma, scale to the my scale. 


The new operator O, at the my scale takes the form 


. 92 zapa 
on = igo mw Sr bg s. (20) 


From the calculations of the corresponding diagrams in Fig. 1 we obtain the Wilson co- 
efficient C'4(my). It is composed of the charged Higgs loop contribution, the chargino 
loop contribution, the neutralino loop contribution and the gluino loop contribution, 


whose analytic expressions are given in the Appendix. 


For the running of the Wilson coefficients including C4 from the my scale to the 
ma, scale we work out the anomalous dimensions by calculating one-loop diagrams 
with operator insertions. We find that all the Wilson coefficients in Eq.(19) run in 
the same way as in the MSSM La, while the new coefficient C4 is not changed, i.e., 
Ca(ma,) = Ca(mw). 


When integrating out A, at the ma, scale we find it gives a contribution ACQ, (ma) 


to the operator Qə 


2 Calma). (21) 


Finally, for the running of the Wilson coefficients from the ma, scale to the m, scale 


the anomalous dimensions are the same as in the MSSM |12]. 


(iii) Case C: Super light A; with mass ma, < mp. 


In this case we retain A; as an active field in the entire analysis. At the my scale 
we integrate out all heavy particles and obtain the Wilson coefficients including Ca. 
Then we run the coefficients down to the mj scale. At the m, scale the effects of O4 
are represented by a change in Co, 


Ó v mem 
AC, = —-—— C 22 


where p is the momentum transfer and I4, is the total width of Aj. Since A; can be 


on-shell in this case, the effects of A, can be sizable even if without tan G enhancement. 


Note that the chargino loop contributions to C4 were also calculated in , where the corre- 
sponding diagrams induced by the A,-squark-squark vertex are neglected since the author 
considered the large tan @ limit. In our numerical calculations we used the full results by 
keeping all terms and thus we also included the diagrams induced by the A¡-squark-squark 
coupling although they contain no leading tan 9 terms. Except for the case of a large tan 6, 
such diagrams induced by the A¡-squark-squark coupling should be included since the 41- 


squark-squark coupling can arise from the F-term of the superpotential and not suppressed 


by the singleness of A;. We checked that in the large tan limit we can reproduce the 


analytical result given in , for the chargino-loop contributions. 


IV. DILEPTONIC B-MESON DECAYS IN NMSSM 


With the effective Hamiltonian and the running of the Wilson coefficients presented in the 
preceding section we calculate the inclusive decays B > X,¢*€~ and their forward-backward 
(FB) asymmetry, as well as the exclusive decays B, — £*(7y. The formulas in terms of the 
Wilson coefficients can be found in |8, |12]. 

Note that our supersymmetric contributions to the Wilson coefficients are given at one- 
loop level (next-to-leading order), while the SM contributions are known at two-loop level 
(next-to-next-to-leading ens In our numerical calculations we consider the one-loop 
results for the NMSSM, while for the SM we also include the two-loop results. 

For the inclusive decays B — X,t we exclude the resonances J/W and V' contributions 
by using the same cuts as in the experiments , ie. the invariant dilepton mass in the 


ranges 
(2m, 2.75 GeV) @ (3.3 GeV, 3.39 GeV) O (3.84 GeV, mp), (23) 


so that our results can be compared with the experimental measurements. 
m and consider the photon in B, — 
(* U^ as a hard photon by imposing a cut on the photon energy E}, which means that the 


For exclusive decays B, — l6 y we follow 


radiated photon can be detected in the experiments. This cut requires E, > 6 mp, /2 with 
ô = 0.02. (Note that for a soft photon both processes B, — 44y and B, — (*(^ must 
be considered together and in this case the infrared singular terms in B, — ¢*¢~y are be 
cancelled by the O(@em) virtual corrections in B, — £*£- [10].) 
In our numerical calculations we perform a scan over the NMSSM parameter space 

2 < tan < 30, — 500 GeV X pers < 500 GeV, 

MER quet 

—50 GeV < Ay € 50 GeV, — 10 GeV < A, < 10 GeV (24) 
with fixed parameters for the sfermion and gaugino sector (500 GeV for all sfermions and 


the gluino, and 200 GeV and 100 GeV for SU(2) and U(1) gaugino masses M» and Mi, 


respectively). In our scan we consider the following constraints: 
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(1) The LEP2 constraints by using the package NMHDECAY lid. 


(2) The constraints from B — X,y which stringently constrain the effective coefficient 


CSF. For the experimental result we use the world average value [20 


Br B= XV = (3.55 £02475, 50:03) x 1074. (25) 


(3) The constraints from B, — " which constrain the Wilson coefficient C4 of the 
"m 


light pseudoscalar operator The experimental result is given by |21 


Br(B, => pty) «1.5 x 107 (90% C.L.). (26) 

Among the relevant dileptonic decays experiment data is only available for Br(B — 
X;u* pu), which is given by [21 

Br(B > Xu u) = (4.3 x 12) x 10$. (27) 


In displaying our numerical results we will show this bound and use it to constrain the 
parameter space. For other dileptonic decay branching ratios, with no experiment data 


available, we will compare the NMSSM predictions with the SM values given by 


Br(B — XTT) = 4.43 x 107, (28) 
Br(B, > yu* u^) = 1.33 x 1075, (29) 
Br(B, => yrr) = 1.85 x 1078. (30) 


Note that the SM prediction for Br(B — X,7+7T7~) was also given in . But our result is 
different from theirs because it is very sensitive to the cuts around the resonances J/w and 

'. While our cuts are chosen as in Eq.(23), we cannot find the corresponding cuts used in 
. We can easily reproduce the result in by varying the cuts. 

In Fig[2] we show the scatter plots of the branching ratio for B — X,p*p” versus tan 6. 
Here we present the results for the three cases: a super light A; (ma, «5GeV), an interme- 
diately heavy A; (5 GeV< ma, < 40 GeV) and a heavy A4 at my scale. In order to see 
how stringent the b — sy constraints are we display the scatter plots with and without the 
b — s^ constraints. From this figure we make the following observations: (1) The branching 
ratio can be greatly enhanced by large tan 8. (2) b — sy constraints are quite stringent 


and can exclude a large part of the parameter space, typically with large tan 8. (3) In the 
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7 6 upper bound y ^ 26 upper bound y ^^ o upper bound | 
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FIG. 2: Scatter plots for the branching ratio of B — Xsu* pu” versus tan 8: the left panel is for a 
super light Aj (ma, « 5GeV), the middle panel is for an intermediately heavy A1 (5 GeV < ma, < 
40 GeV) and the right panel is for a heavy A4 at my scale. The dark (red) points are allowed by 


b — sy, while the light (sky-blue) points are excluded by b — s». 


parameter space allowed by b — sy the decay B — X¿u* pu” can still be greatly enhanced, 
especially in case of a super light A. The 20 experimental bound on B — X;u*gu- can 
further exclude a large part of the parameter space. Almost no points with tan 8 > 15 in 
the parameter space survive all the constraints. From the left panel of Fig[2] we see that 
some part of the parameter space with a super light A, is still allowed by b — sy and 
Beg uu. 

Let us take a look on the constraints from the process B, — ~*~, whose branching ratio 
is given by 
2 


Vis 
0.04 


2 3 
Br(By = tu) = 1.2.x 10-7 | 72 | fm == 
i D ü | mo 


1.49ps | |245MeV 


2 
x [C3, + (Co. " 2 C9 (31) 


We see that the contributions are from Cio and Cg,,. While the contribution from Cio is 
suppressed by the factor m,/mp,, the contributions from Cg,, can be enhanced by large 
tan 3 (Co, , contain terms which are proportional to tan? 8, as shown in the Appendix). 
This feature can be seen from Figl3] in which we set aside the b — sy constraints and 
illustrated the B, — utu constraints. We see that, similar to the b — sy constraints, the 


B,  u*p- constraints are stringent for a large tan 8. If we impose the b — sy constraints 
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which exclude a very large tan 5, then the further constraints from B, — ^p^ are stringent 


only for the parameter space with a very light Aj. 


2c upper bound | 


26 upper bound | 26 upperbound | 


Br(B> X,u* w) 
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tan tanp tanp 


FIG. 3: Same as Fig[2] but with the B, — tp constraints. The dark (red) points are allowed 


by B, —^ wtp, while the light (sky-blue) points are excluded by B, > pty". 


'The results for other dileptonic decay branching ratios, for which no experiment data are 
available, are presented in Fig] To see how stringent the constraints from B — X,u*u- 
are we display the scatter plots with and without such constraints (all the points satisfy 
b — sy and B, > pty"). 

From Figli] we see that under the constraint from B — Xyu*u”, the branching ratio of 
B — X,r*r- does not deviate significantly from the SM value. The reason is that these 
two decays are highly correlated except for the contributions of C, , and C4 which are 
dependent on the lepton mass. If the contributions of Co, , and C4 are dominant, then 
Br(B — X,r*r ) should not be severely constrained by B > X,utp~. As discussed below 
Eq.(31), the contributions of CQ, , and C4 are important only for very large tan 9 which is 
not allowed by b — sy. As a result, the contributions of Co, , and Ca are not dominant and 
thus B > X;,r*T- is highly correlated to B > X u” u”. 

Note that for B > X;,r*r- it may be rather challenging to disentangle the NMSSM 
effects from the SM value in future experiments. One reason is that, as discussed above, the 
NMSSM effects are no longer so sizable under the constraint of B — Xsu*u”. The other 
reason is that the SM prediction has its own uncertainty. If we consider the uncertainty of 


the input SM parameters, we can obtain the uncertainty (about 2096 as found in ) of the 
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FIG. 4: Same as Fig[2] but for B —> X,r'*r , B, —> yutp and B, => yr*vr- 


(red) points are allowed by B => Xspu*p”, while the light (sky-blue) points are 
Bo Xsu™ u. 


The dark 


excluded by 


SM prediction. But in Fig] we did not show such an uncertainty of the SM value because 


for all the results, both NMSSM and SM, we used a same set of the SM parameters without 


allowing them to vary in the uncertainty range. Since the SM parameters are involved 


in both the NMSSM and SM values, all the results are subject to some uncertainty if we 


consider the uncertainty of the SM parameters. Of course, such uncertainties will deteriorate 
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the observability of the NMSSM effects. 


— m, 248.889 
L 1 ] 
0.6 | w 4 


FIG. 5: The forward-backward asymmetry of B > X,u* p^ versus s = p?/m? (p? is the invariant 


mass of u* and 7). 


Finally, in Fig[5] we show the results for the forward-backward asymmetry in B — 
X,u*pn- under the constraints from B —> X,u*u-. 


V. SUMMARY 


In the framework of the NMSSM we examined the rare dileptonic decays B — X,0*£^ 
and B, — (*(7 y paying particular attention to the light CP-odd Higgs boson. We found 
that in the parameter space allowed by current experiments, such as LEP II and b — sy, 
the branching ratios of these rare decays can be greatly enhanced and thus the experimental 
data on B > X,u*p further stringently constrains the parameter space, especially with 
a super light CP-odd Higgs boson and large tang. In the surviving parameter space we 
gave the NMSSM predictions for other unmeasured dileptonic decays which may hopefully 
be measured at the future LHCb or super B factory. 
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APPENDIX A: WILSON COEFFICIENTS 


The Wilson coefficients Cz, Cy and Cio in the NMSSM are the same as in the MSSM 
. Here we present the new coefficient C4, and Cg, and Cg,, whose predictions in the 

NMSSM are different from the MSSM. We checked that we can analytically reproduce the 
MSSM results El n" (However, the NMSSM results can not explicitly reduce to the 
MSSM results by simply dropping out the singlet $ (say setting À = k = 0) because the 
y-term is generated by $). 

Although in our numerical calculations we used the complete results by keeping all terms, 
here, for simplicity, we only present the terms which can be enhanced by large tan 8. At the 
my scale each Wilson coefficient is composed of the charged Higgs loop contribution from 
Fig.1(a-c), the chargino loop contribution from Fig.1(d-g), the neutralino and gluino loop 
contribution from Fig.1(h-k). 


For the charged Higgs contributions: 
¡MAA 


a =— tan OF) (xg, xwt) ay) 
goamw 
2 2,42 
E mym MH ae 
Cg, =- iy tan? g| EU UP gie ewe) + 33 Fi (tigt, ew) eo 
Am miy Mwy mM yH 
2 
+ mpm mM Ht 
CE = tan? b|| 00 bet” | Ales £w) 
Aa Ww 
mem, 
+ z Pi (Ley 5 iw) (A3) 
Ty TE E 


Note that although H*W+A, vertex has 1/tan/ suppression by the singleness of Aj, 
H*G* A, vertex which comes from the soft term AAS HIH? does not have such 1/ tan 8 


suppression. Thus the contribution to C4 from the loop involving H* and G^ with H*G* A, 


coupling (in the Feynman gauge) is proportional to tan 6. 
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For the chargino contributions: 


ye _ „tang toe 1/ 
Ci = iS Goad) Ó_ bu x vow? (2 hax? c) 
1/2 * 
= Lm C — Rui Fo Ga ixi sese) | H (A4) 
VU. Tuc 
a+ Lm: al 
CS = an? 8 Y > Di (,k,3,1 A tit 
i,k=1 j,l=1 
EET 
Sy Oir ti1Xf> »Tzoze) + (v; iiXf? psg? 
M 
b — bg ti DEL 
tk—1 
mj 2j 
+ Ok Z” ZU Fs, NF E  BCxC, £ 0x0) 
Xy 
TROC ou. qu 
_ Xr ZA s Fs(zj, af TEZE, By xc) H (A5) 
myc 
3 2 V2U4 UA m Be 
xt Mam] 2 ; a al’ al 
CX = — tan”? B I'4(i, k, 7,1 A 0 a 2 dz LEO 
34 Am, 22 a pm MED 
2204 Mye 
LESE |- 2i; Fo (Eg, ET. sese) + Raji Fig, ET I2 
92 Wie 
V2UA T? mm, c 
J 
ee EU p TCin) 
tk—1 
2 
MA, 2j 1l 
t—3204 Z7 ZY Firg 13G, Eege Toge) 
X 
MIO qu 
——XA ga Ze F3(2;, axe TREO To zl (A6) 
Tac 


Note that here Cg,, contain terms which can be enhanced by tan? (the overall factor 
tan? 8 multiplied by 1/ cos 8 gives tan? 8 in large tan limit). 
For neutralino contributions: 


0 i  tanf * * m 
CX =- NUM phun [gne oF 
A VgVa 008 Ow 1 ue bg fun 


i— 1x3 ? vios) 


Be Fox. x9 Tx 0x 2) a l2(i, j, Era, Tan, 7 


ó. cot Ü v 1/2 
o, .) + dik? === -I s x ol alts j, EX CN } (A7) 


xF (v; 


i—1bk'—1 
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Es 2v. B V2UH UH m. o 

y 1 mpm Y Y Y al^ all 
X = eee RM f t n? 8 N! a EUN. iT ki P T7 ~ 
Y l KyKi, Am, cos Oy 1 my wT a(t, J, E) PL hax) 


20E ZIK | aah 
n lori D Fox 1x9 T) 


+ Qa (= e DEE, , m 
Mya xj 
Xj kl 773j* rpil* rpk'2 
301 AT Tai iK) + T ZY TE TRO FG, iby o Egia) 
bi 
2 Tn. o 
m X : x * 35* ee 
e ( ? «(i k) zU 9m Zim gak ys EL jx Bx 050, iso) 
mo LM 
(razo mz "EM E Pd), q ttt 7: (A8) 
a 5 6 VIUA Ud m 
EU 1 Mamy 2 aU. a 
QE Sia AE Ix A e S PERS Sa 
9.7 akp mmh, cosy OP De de UN |g eT IPG) 
2UA Z3F [q 
Be [Rarer PU 1x97 75070) 
ERES. C— XE TIS "p AE SrH ons m 
E " 


My 
A ak! sum op £u - 
+204 T3 To(i, j,k ICAA 


bki 
m? my * * 35" v 
Ec G j Te(i, zm ZEE zs O) Fa(vpo 1599 ERIK x x x9 0) 
x9 LN 
k 
(razo mz Dies «(Fig digs 2h (A9) 


For gluino contributions: 


puc LU 
n 6 "gb 


lji 
+T3 ji (Lo, a Tw, |) ZR RI (A10) 


5 i 89% 1 
C% = 93 tan BTS D Pi (z; bi ig) 


zt 
Veo Ves 3g 92 


z 1 mpm T 2 | U AD i 
Ch. = TST z oj; Pi (xz 
Q^ KpKi 4 3gim aes my ^7 1r, ag) 
Ue 
te Fi (xz, by > Ta, +) , (A11) 
bj-1 
3 
a 1 Be 16g2m; 2 8 rmjl* DS Uâ ES 
(UL. = — 32,2 tan TËT? ôi; Pilaz = 
UA pos 
CRUS Fs RR as (A12) 
b5-1 
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In the above expressions the constants and functions are defined by 


du: ste T 
N; = ZZN sin 0, — Ze cos Oy, Nj —— y sin 0, + ZA cos Oy, 
A 2j 
Rotj = — (UA ZIZ -UAZIZE) UAI 
~~ B v2 ° 2: 
A 1 z2j 
Q "n 22 D H2 7 a +UE Hyg 2 e E H gaz H 
"T s m 
HH ES 1 44 75 34 
HEC -Rir = = a Vlast, (ZN EN; tZ MN) + V2X(Z MZR + ZN ZN) 
$US [> (ZEN; + ZEN) — VOM ZN ZN + ZA ZE) 


"TE 
—V2kUA. pcr qu ZAR) 


QU, = QR, = su Ace INN; + ZNN) + VIM ZN ZN + ZN ZN) 
EN " FN; + ZR NL) + V2X(ZN ZN + ZN ZR] 


—V2kUE (ZŠ ZRF + Z ZRF) 


Ti(i,k, 5,1) = (TPT — bid) Z!" Z2 — S ui ea D 
V2my sin " 
To(i, j,k’) = V2tan 6, 2 TU TE3 ZA pH TE? 
v2mw cos 8 
T Et Th T} 
sli, j) puro A 
Poli) NTE TBE ZI" (ZW sin Ou + Zi cos 0u) ZF 
cos 
cnp Um quA rar OS] 
Te = A3 go gam gem = ATE TET ETE UE 
mem + Ti Ti) + (ATË Te? + ATTE TEP) 
+ A5 (T TT + TÈT, TED) + AT? TES 


Tm = —2m,UH TER 20e TE ps mc + (ATS T + ALTE UM 


TIP = ATTE” — As(TBTE. + TET.) 


À A 
A] = yg Us is) ADU, Æ = gos + sU) + ADUE 

A A 
As = 77 (va Uds + Ucn) — Avda, Aa = To(va as + SUIT) + AvUad 

z 4 , 2 
= seg E 3 Sin? 0w YU Ag = q tan? OS 
A; = mpeg tan? QUE. Ag = (1 + 2eq sin 0 US 
2 SYM 
zlnz 1 rlnix ylny 
P. = —— F = — 
(2 295. Alay = (29-19) 
1 r?lng y?lny 
F: = — 
I 

F3(x, y, ) — see 


7 7 (r- D - yy - z) 


22:2 vi. 
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